% График напряжонности поля на тарелке
Angls=360;
for n=-n_1:n_1
  ro(n+n_1+1)=R0/n_1*n;
end
clg;
plot(ro(1:2*n_1+1),Frup(Angls/2-n_1+1:Angls/2+n_1+1),'m-');title('эпюра напряжонности');
  xlabel('Радиус тарелки,в см.')
  ylabel('напряжонность поля')
axis([-R0 R0 0 1]);
pause
close
disp('Отношение напряжонности на краю зеркала к напряжонности в центре равно')
f_1/Frup(Angls/2+1)
pause

